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ABSTRACT 

Context. The synchrotron X-rays can be a useful tool to investigate the electron acceleration at young supernova remnants (SNRs). 
Aims. At present, since the magnetic field configuration around the shocks of SNRs is uncertain, it is not clear whether the electron 
acceleration is limited by SNR age, synchrotron cooling, or even escape from the acceleration region. We study if the acceleration 
mechanism can be constrained by the cutoff shape of the electron spectrum around the maximum energy. 

Methods. We derive analytical formulae of the cutoff shape in each case where the maximum electron energy is determined by SNR 
age, synchrotron cooling and escape from the shock. They are related to the energy dependence of the electron diffusion coefficient. 
Next, we discuss whether information on the cutoff shape is provided by near future observations which gives simply the photon 
indices and the flux ratios in the soft and hard X-ray bands. 

Results. If the power-law index of the electron spectrum is independently determined by other observations, then we can constrain 
the cutoff shape by comparing theoretical predictions of the photon indices and/or the flux ratios with observed data which will be 
measured by NuStar and/or Astro-H. Such study is helpful in understanding the acceleration mechanism. In particular, it will bring us 
another independent constraint on the magnetic field strength around the shocks of SNRs. 

Key words. Acceleration of particles - ISM: cosmic rays - ISM: supernova remnants - Turbulence - X-rays: ISM 
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1. Introduction 

Observations of synchrotron X-rays tell us that the young su- 
pernova remnants (SN Rs) are the ele ctron accelerator up to 
10-100 TeV (e.g.. lKovama et alJll995h . Spectral fitting of the 
synchrotron radiation from radio to X-ray bands gives us in- 
formation on the acceleration mechanism. For example, the 
maximum energy of electrons £ maiie and the magnetic field 
strength B are constrained by measured roll-off fre quency of 
the spectrum, v ro u, which is propo rtional to BE m . dx 2 ( Re vnoldi 
1998: [Reynolds & KeohaneH l999). Furthermore, if rapid cool- 
ing time of 1-10 yr is responsible for the observed narrow ri m 
brightening dVink & Laming||2003l: iBamba et al1l2003l l2005allbT: 
Yamazaki et al. 2004) an d/or time variability of the synchrotro n 
X-ravs (lUchivama et all 120071: lUchivama & Aharonianl 120081) . 



show that observations by these satellites will play an important 
role in this context ( § [3] and |4j . 



2. Cutoff shape of electron spectrum 

In this paper, we assume the shape of the electron spectrum 
around the maximum energy, E max<e , in the form, 



N(E)^E-"exp[-(E/E m ^ e ) a ] 



(1) 



then B ~ 0.1-1 mG is inferred, so that E„ 



10 TeV. Given 



the values of E, mx e and B as well as the SNR age f ; 
even estimate the maximum energy of protons E maXtP 
proportional to E mlix , e 2 B 2 t Rge (see § IA. lb . 



age, we can 
which is 



F or young SNRs, Vmii is ty pically below the X-ray band 
(e.g., Bamba et al. 2003, 2005a b), so that synchrotron X-rays 



are emitted by electrons whose energy is near E miXfi . Hence it 
is expected that information of the electron spectrum at high- 
est energy range can be extracted from the soft and hard X-ray 
data. In particular, recently launche d satellite The Nucl ear Spec- 
troscopic Telescope Arr ay (NuStar) (Hailev et al.ll2010 ) and near 
future mission Astro-H (Takahash i et al.ll2010h will observe hard 
X-rays whose photon energy is larger than 10 keV. In this pa- 
per, we focus on the cutoff shape of the electron spectrum and 
show that it may provide another independent way to tackle the 
problems of particle acceleration at young SNRs (§|2). Then, we 



As discussed in the following, the spectral index p and the cutoff 
shape parameter a contain rich information on the acceleration 
mechanism. 



2.1. Spectral index p 

When the diffusive shock acceleration works at the adiabatic 
shock with a compression ratio of r and the energy loss effects 
are negligible, the spectr al index is given by p = (r + 2)/(r - 1 ) 
in the test-particle limit dBelllll978: Blandfo rd & Ostriked l978h . 
In particular, p = 2.0 in the case of the strong shock limit 
r = 4. This simplest case may not be for actual young SNRs, in 
which GeV-to-TeV gamma-ray observations of youngest SNRs 
such as Cas A and T ycho suggest p > 2.0 dAbdo et al.ll2.Q10b 
iGiordano et al.ll20T2h . This fact has been already inferred both 
from radio synchrotron spectrum and fr om propagation model 
of cosmic rays from SNR to Earth (e.g.. | Strong & Moskalenkol 
ll998blPutze et al.ll2009tlShibata et allfeoi ll) . 

If the magnetic field is strong enough (e.g., in the Bohm dif- 
fusion model, K(E) cc E, see equation dA.51 >). the synchrotron 
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cooling is responsible for a spectral break, above which the elec- 
tron spectrum becomes so fter and the sp ectral index increases by 
1.0, so that p > 3.0 (e.g.. lLongairll 1 9941) . Note that sometimes 
the synchrotron cooling effect causes spectral hardening or pile- 
up dLongarrill994t iDrurv et al.|[l999l IZirakashvili & Aharonianl 
l200lT For typical parameters of young SNRs, however, such 
hard component can hardly be seen for more than a decade of 
electron energy. At most a small bump is formed just below 
£max,e, which may be seen as a small excess in the radiation 
spectrum (see §[4]i. 



2.2. Cutoff shape parameter a 

The cutoff shape parameter a also depends on the details of elec- 
tron acceleration, such as the magnetic field strength and the en- 
ergy dependence of the diffusion coefficient. In this paper, we 
assume that the diffusion coefficient of high-energy electrons in 
the following power-law form: 



K(E) oc E 13 



(2) 



Usually the Bohm diffusion, ft — 1, is widely adopted. How- 
ever, in general ft may deviate from unity. For example, it is 
well known that ft becomes 1/3 if the particle diffusion is con- 
sidered in the Kolmogorov ma gnetic turbulence, while ft - 1/2 
for Kraichnan turbulence (e.g., Blandford & Eichler| [l987Pl So 
in these cases, the value of ft even tells us the properties of 
magnetic turbulence. In another context, if the wave damping 
due to ion-neutral collisions is sig nificant, ft may approach m 2 
dBvkov et alJliOOOllLee et ai1l2012h . 

In the cooling limited case, where E m - dXfi is determined by 
the balance of acceleration and synchrotron cooling (t dcc (E) = 
t syn (E)), the spectral shape factor a is related to ft as 



■ ft+l 



(3) 



which is analytically derived as in § IB.1I In the age-limited case, 
where synchrotron coolin g effect is neglect ed and E maxfi is de- 
termined by the finite age. lKang et alJ d2009t) have given 



a = 2ft 



(4) 



by fitting their result s of numerical simulation (see also 
iKato & Takaharal 120031) . Note that in the case of Bohm diffu- 
sion (that is, ft — 1), both age-limited and cooling-limited cases 
give the same value, a — 2. If ft + 1, the values of a for the two 
cases are different. 

It may hap pen that the maximum energy is limited by the 
escape process fPtuskin & ZirakashviTill2005t IDrurv et ail [2009: 
Capr ioli et al.l 12009b lOhira et alJl201Qp . Here we consider the 
simplest case (see § IB.2t . in which a free escape boundary ex- 
ists upstream of the shock front. In the test particle limit, we 
analytically derive 



a=ft 



(5) 



If nonlinear effects, in particular the decay of self-excited up- 
stream turbule nce, is taken in to account, then a may be slightly 
larger (e.g., iLee et al.l l2012h although precise cutoff shape is 
at present highly uncertain. However, as long as ft m 1, the 



In general, if the spectrum of magnetic turbulence is in a form of 
E k oc k~ s , then ft and i are related as ft = 2 - s under the assumption that 
accelerated particles scatter via wave-particle resonance interaction. In 
particular, for Kolmogorov (s = 5/3) and Kraichnan (s = 3/2) turbu- 
lence, we obtain ft = 1/3 and 1/2, respectively. 



escape process doe s not affect the maximum electron energy 
(Ohira et al. 2012b) and hence the cutoff shape for young SNRs. 

So far, we have discussed the cutoff shape in the diffusion ap- 
proximation for particle motion. Around the maximum energy, 
a particle mean free path tends to be longer so th at the particle 
motion becomes ballistic (e.g., IVladimirovll2009l) . Our discus- 
sion could be extended to the case of superdiffusive transport, 
whose o bservational evidences ha ve been found in interplanetary 
shocks dPerri & Zim bardo 2012|). In such cases, cutoff shape 
may be modified. It is expected that the highest energy parti- 
cles go out of the acceleration region more rapidly than the case 
of diffusive transport, so that the spectral slope may be steeper. 
Anyway, it can be said that cutoff shape parameter a contains 
rich information on the particle acceleration mechanism. 



3. Synchrotron X-ray diagnostics of cutoff shape 

In this section, we consider the synchrotron radiation from elec- 
trons whose energy distribution is given by equation dT}. As- 
suming an isotropic pitch-angle distribution of electrons, the 
energy spectrum F v [erg s _1 cm _2 Hz _1 ] of the synchrotron ra- 
diation is computed. The energy flux in the photon energy 
range between ei = hv\ and ei = hv2 (vi < V2) is given by 

XV? 
" F v dv . We calculate photon indices and flux 

ratios in various energy bands, on which we focus in the fol- 
lowing of this paper. These spectral quantities are determined 
if we specify four parameters a, p, E mdXfi , and the magnetic 
field strength B in the emitting region. Note that the flux nor- 
malization is not necessary. The field strength B is only re- 
quired to determine the frequency that gives the peak of the 
vF v of the synchrotron radiation, so called a roll-off frequency 
VroiiO* BE mdXy 2 ), which roughly corresponds to the character- 
istic synchrotron frequency of ele ctrons with E m . dx , e dRevnoldsl 
1998; Re ynolds & Keohandf 19991) . Different parameter sets but 
with the same value of BE mdx 2 give the same photon indices 
and the flux ratios. Therefore, independent parameters are a, p, 
and BE mdx /. 

The left panel of Figure [TJ shows 10-50 keV photon index as 
a function of 2-10 keV photon index. We adopt p = 2.0 (thick 
lines) and 3.0 (thin lines), and a = 0.5 (light blue), 1 (green), 2 
(red), 3 (blue) and 4 (purple). Along each line, both a and p are 
constant, and the quantity BE mdXfi 2 changes. In the limiting case 
of BE mdXy 2 — > 00, rolloff frequency v ro n goes beyond the ob- 
servation band, which implies the synchrotron spectrum is well 
approximated by F v oc y-b'-V/ 2 . Then the photon index becomes 
an asymptotic value, 



p+l 



(6) 



so that F M = 1.5 and 2.0 for p = 2.0 and 3.0, respectively, which 
correspond to the left end of each line. When BE mdx , 2 decreases, 
the rolloff frequency v m \\ crosses the observation bands, so that 
the photon indices become larger than the asymptotic value r^. 
After passing through the harder band 10-50 keV, v ro u crosses 
the softer 2-10 keV band, so that the photon index in the former 
band is larger than the latter. For fixed p, if a becomes larger, the 
flux beyond v m \\ more rapidly decreases, resulting larger photon 
index. Hence each line has steeper slope for larger a. In this 
paper, the value of BE mdx e 2 is reduced to !Q 3 pG (TeV) 2 , which 
corresponds to the characteristic frequency of the synchrotron 
radiation of 2.9 x 10 16 Hz. Light blue lines in all the Figures as 
well as the thick green line of the right panel of Figure |2]have 
the right end, which corresponds to this lower limit. 
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Fig. 1. Photon index in 10-50 keV band (left panel) and the flux ratio R, = F(10 - 30 keV)/F(30 - 80 keV) (right panel) as functions of 
2-10 keV photon index. Along each line, both a and p are constant, and the quantity BE mlixs 2 changes from l0 3 pG (TeV) 2 to <x>. Thick lines are 
for p = 2.0, while thin lines are for p = 3.0. Light blue, green, red, blue and purple lines are for a = 0.5, 1, 2, 3 and 4, respectively. 



The right panel of Figure[T]shows the flux ratio R\ = F(10 - 
30 keV)/F(30 - 80 keV) as a function of the 2-10 keV photon 
index. In the case of BE m . dXte 2 — > oo, the hardness ratio R = 
F(si - S2)/F(s3 - £4) of two energy bands s\ - £2 [keV] and 
£3 - £4 [keV] has an asymptotic value: 



2-r„ 



2-r„ 



R 



2-r« 



2-r„ 



(7) 



for Too + 2, while R — > ln(£2/£i)/ ln(£4/£3) for F M = 2. Hence 
in the present case (e\ = 10 keV, £2 = £3 = 30 keV, and £4 
SO keV), we have r«, = 1.5 and /?, -> 0.668 for p = 2.0, while 
Too =2.0 and R\ — > 1.12 for p = 3.0. When BE m!iXfi 2 becomes 
small, both the flux ratio and the 2-10 keV photon index become 
large, however the decay slope is steeper for large a. 

The left panel of Figure [2] shows the hardness ratio R2 = 
F(2 - 10 keV)/F(10 - 80 keV) as a function of the 2-10 keV 
photon index. In the limit of BE m . dx ^ 2 — > 00, we have r M = 1 .5 
and R 2 -» 0.302 for p = 2.0, while r m = 2.0 and fl 2 -» 0.774 
for/? = 3.0. 

The right panel of Figure [2] shows the flux ratio Ri = 
F(10 - 30 keV)/F(30 - 80 keV) as a function of R 3 = F(2 - 



10 keV)/F(10 - 30keV). In the limit of BE n 



00, we have 



Ri -> 0.668 and R 3 -> 0.755 for p = 2.0, while R { -> 1.12 and 
R 3 1.46 for p = 3.0. 

Both panels of Figure [3] are the same as Fig. Q] but for 
p = 2.3 and 3.3. The former is typical for the source spec- 
trum of Galactic cosmic rays, which is inferred by the propaga- 
tion model (e.g..lStrong & Moskalenkolll998UPutze et alJl2009t 
IShibata et al.ll2.01lh . It is also expected from gamma-ray obser- 
vations that young SNRs such as Cas A and Tycho have th e en- 
ergy spectrum with p = 2.3 (see § I4T21 1 Abdo et alJ feOlO). The 
latter case (p - 3.3) is realized if the synchrotron cooling is 
significant enough to provide a cooling break below which the 
electron spectrum has p = 2.3. In the limit of BE mix e 2 — > 00 , we 
have r M = 1.65 and R { -> 0.779 for p = 2.3, while r M = 2.15 
and Ri — > 1.31 for p - 3.3. Since differences are too small, we 
do not show the counterparts of Fig.|2]for p = 2.3 and 3.3. 



4. Discussion 

The value of p may be determined by the spectral slope of the ra- 
dio synchrotron and/or GeV-to-TeV gamma-ray spectrum. Once 
p is independently determined, then comparing theoretical lines 
of Fig. [T]or[3] with observed data which will be provided by NuS- 
tar and/or Astro-H, we can constrain the value of a. As discussed 
in § 12.21 it is helpful in studying if the acceleration is age-limited, 
cooling-limited or escape-limited. In particular, we can check if 
a 2 or not, which corresponds to the popular case of Bohm 
diffusion in age- or cooling-limited acceleration. One may ex- 
pect that a = 2 is typical. Indeed, some observational results ar e 
consistent with the case of a - 2 (e.g., iNakamura et al.ll2012l) . 
although no firm conclusion has been derived. 

In addition to other observational information such as wide- 
band radiation spectrum, X-ray spatial and temporal variability 
and so on, the value of a will bring us another independent con- 
straint on the magnetic field strength (for typical example, see 
§ 14.1. 11 1. In another point of view, if we know, by other observa- 
tions, the field strength and know how £ max . e is determined, then 
with the aid of Eq. OJ, or (0, we can further constrain the 
value of B, which brings us rich information on the acceleration 
process, in particular, the electron transport mechanism. 

Note that our present method is able to constrain the value 
of a without detailed spectral fitting which needs brightness. It 
is easier than the spectral fitting analysis to observationally de- 
termine the flux ratios and/or photon index. Even for sources 
whose X-ray brightness is too small to perform the precise spec- 
tral fitting, we will be able to discuss the spectral parameter a. 
Hence we can obtain larger sample, which enables us to do the 
statistical discussion. Consequently, we will be able to extract 
more general properties of particle acceleration at SNRs, which 
cannot be done by individual analysis for the small number of 
bright sources. 

The left panel of Fig. [2] shows that all the curves are de- 
generated with each other, along which the observed data points 
should lie if the X-ray emission arises from the synchrotron ra- 
diation. If there are outliers, we can discuss the existence of 
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Fig. 2. The hardness ratio R 2 = F(2 - 10 keV)/F(10 - 80 keV) as a function of the 2-10 keV photon index (left panel), and the hardness ratio 
R, = F(10 - 30 keV)/F(30 - 80 keV) as a function of R 3 = F(2 - 10 keV)/F(10 - 30 keV) (right panel). The meaning of each line is the same as 
Figure 1. 
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Fig. 3. The same as Fig.[T]but for p = 2.3 and 3.3. 



extra components, such as a bump of the electron spect rum near 
£ ma v e due to effe cts of pile-up via synchrotron cooling (jLo ngair 
1 994; iDrurvet alfl999t IZirakashvili & Aharonianl2()07b and/or 
nonlinear acceleration (e.g., Malkov & Drurv 2001) , jitter or 
diffus ive synchro tr on emission (j Toptygin & Fleishman 1987; 
Med yedevl 12005 iReville & Kirk I 1201(3: iTeraki & Takaharal 
20111) . secondary synchrotron radiation origin ated in acceler- 
ated protons generating charged pions (e.g.. lYamazaki et al.l 
2006 ), nonthermal bremsstrahlung emission ([La ming 20 01a1fbt 
Vink & Lamingl2003kl^]d2008t(ohira et dJ2012«£ and so on. 
Similar discussion may be done with the right panel of Fig. [2] 

A caveat is that one-zone model is adopted for simplicity. 
Clearly our plots change if we see emissions from different re- 
gions which have different parameter set of B, E m . dXfi , p, a as 
well as the flux normalization. In order to avoid or minimize 
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this problem, hard X-ray observations of small emission regions 
are necessary. In this sense, hard X-ray imaging system onboard 
NuStar or Astro-H is useful. 

In this paper, we have focused on hard X-rays which arise 
from synchrotron radiation of electrons near E m ^ e . It may be 
possible to do similar analysis in gamma-ra y bands using nex t 
generation gamma-ray telescopes like CTA (lActis et al.l 1201 ll) . 
However, at present there are several uncertainties. First, the 
emission mechanism is uncertain: even for representative case 
of SNR RX J17 13.7-3946, we have not yet determined whethe r 
the gamma-ray emission is leptonic or hadronic (see § 14.1.11 . 
Second, in the leptonic case, where the gamma-rays arise from 
inverse Compton emission of accelerated electrons, seed pho- 
tons are uncertain; in addition to CMB, infrared and optical pho- 
tons, whose intensity is uncertain and depends on the position 



R. Yamazaki et al.: Electron cutoff shape 



of SNRs in the Galactic plane, significantly contribute to the 
gamma-ray emission. Finally, systematic error of the gamma- 
ray measurement might be too large for our discussion. If all of 
above uncertainties go out, the flux ratio and/or photon index in 
gamma-ray bands will be another independent diagnostics. In 
particular, if gamma-rays are hadronic, we can obtain the cutoff 
shape of protons near E max p , which cannot be done by hard X- 
ray analysis proposed in this paper. However, even in this case, 
we need a careful analysis because spatial inhomogen eity of tar- 
get matter modifies the gamma-ray spectrum (e.g., lOhira et all 
l201llllnoue etaTEoTl . 



4.1. Application to specific cases 

In the following, let us consider the case of Cas A and 
RX J1713. 7-3946 as representative examples. 



4.1.1. RX J1 71 3.7-3946 



Tana ka et al.l (|2008) presented results of Suzaku observation of 
RX J1713.7-3946. XIS and HXD onboard Suzaku detected X- 
rays from this SNR in 0.4-12 keV and 12^4-0 keV bands, and 
measured photon indices in these bands are 2.39 + 0.01 and 
3.2 + 0.1, respectively. We adopt these values as 2-10 keV and 
10-50 keV photon indices, re spectively. Furtherm ore, using the 
results of spectral analysis of lTanakaet al. (2008), we calculate 
the energy flux of the whole SNR as F(2 - 10 keV) = (4.59 ± 
0.04) x 10~ 10 erg s-'cirr 2 and F(10 - 80 keV) = (1.32 ± 0.05) x 
10~ 10 erg s _1 cnr 2 , which lead to the flux ratio, R 2 = 3.5 ± 0.2. 
Derived value of R2 and the adopted 2-10 keV photon index lies 
on theoretical lines in the left panel of Fig. [2] which implies that 
the X-ray emission is the synchrotron radiation. 

Recently Fermi measured the gamma-ray spectrum of 
RX J1713.7-3946 in 3-300 GeV energy range, whose photon 
index (1.5 ± 0.1) is typical for leptonic inverse Compton emis- 
sion from high-en ergy electrons with the spectral index p « 2.0 
dAbdo et al.ll201 lb . Hence, at first, let us consider the possibil- 
ity that gamma-rays are emitted by inverse Compton scattering. 
Assuming CM B and infrared photons as seeds of the scattering, 
lLietal.l(l20Tl extracted the electron distribution from the ob- 
served gamma- ray spectrum and obtained p = 2.0 and a = 0.6. 
iLee etal.1 (120121) also obtained similar result, a = 0.5, to explain 
the gamma-ray spectrum. The value a = 0.6 implies B = 0.3 
in the age-limited case, while B = -0.4 in the cooling-limited 
case. The latter is unlikely because negative value of B is quite 
unnatural in the context of the magnetic turbulence. This is 
also inferred from the fact that the flux ratio of synchrotron X- 
rays to inverse Compton gam ma-rays leads to the field strength 
B ~ 10-20 yuG (e.g.. iKatz & Waxmanl l2008t I Yamazaki et alj 
2009; Ellison et al. 2010), so that synchrotron cooling effect is 
not significant- However, one can find from Figs. Q] and |3]that 
no lines for a = 0.5 are consistent with the measured photon in- 
dices in 2-10 keV and 10-50 keV band s. In particul a r, a » 1.5 
if p x 2.0. This fact was indicated by iTanaka et al.l d2008l) . in 
which they found that the X-ray spectrum taken by Suzaku is 
consistent with both a — 1 and a = 2 cases. Therefore, the 
simple leptonic inverse Compton model fails to explain both the 
X-ray and the gamma-ray spectral shapes simultaneously. 



2 There are several possibilities to explain observed time variability 
and narrow rim brightening of the synchrotron X-rays w i thout ampli- 
fication of the magnetic field (e.g., Katz & Waxman 2008; Bvkov et al. 
2008). 



On the other hand, hadronic scenarios for the observed hard 
gamma-ray spectru m are still viab le if we consider the shock- 
cloud interaction dlnoue et al.ll2012l) or extreme limit of nonlin- 
ear particle acceleration (e.g., I Yamazaki et al. 2009). In these 
cases, the strong magnetic field (B ~ 0.1-1 mG) is predicted, 
which is also inferred from the detection of time variability 
and narrow rim brightening of the sync hrotron X-rays although 
there are some counter arguments (e.g. , IKatz & Waxm an 2008; 
iBvkov et alj2"00 8; Reynolds et al.l2012l) . Hence, the electron ac- 
celeration is limited by synchrotron cooling, so that the spectral 
index of electrons near E maxfi should be p > 3.0 (see § 12.11 ). One 
can see from Fig. [T]that if p m 3.0, then the observed photon 
indices in 2-10 keV and 10-50 keV bands tells us a > 2, so that 
B > 1 . Therefore, Bohm diffusion is consistent with observation. 

Note that, however, this argument is not conclusive. We con- 
sider the spectrum of the whole SNR. As discussed previously, 
one-zone model may not be well. Further observations resolving 
smaller emission regions are necessary. 



4.1.2. Cas A 

XIS and HXD onb oard Suzaku have m easured the spectrum in 
3.4^10 keV band dMaeda et al.l 12009). They found that non- 
thermal component in this energy range is described by a sin- 
gle power-law form with a photon index of 3.06 + 0.05. Here 
we adopt this value as 2-10 keV photon index, although a ther- 
mal component prevents us from determining the index pre- 
cisely. Using their fitting parameter, we have F(2 - 10 keV) = 
(6.5 + 0.9) x 10~ 10 erg s"'cm -2 . Other observ ations in higher en- 
ergy b ands such as Swift/BAZT (14-195 ke V) (|Bau mgartne r et al.l 
2012) and BeppoSAX/PDS (15-300 keV) ( Vink & Lamingl2003l) 
show softer photon index of 3.26 + 0.09 and 3.32 ± 0.05, re- 
spectively. Hence we expect the photon index in 10-50 keV 
band is slightly softer (3.1-3.3), than the value 3.06 measured 
by Suzaku. Using the parameters given by 5w//t/BAT, we have 
F(10-80keV) = (9.6±0.7)xl0" n erg s^'cirr 2 '. Then we obtain 
the flux ratio R 2 ~ 6.8 ± 1.6. 

Although the index of the electron spectrum, p, is not fixed 
yet, measured gamma-ray spectrum shows that p is larger than 
2.0 dAbdo et al.ll2010b . Here we adopt p « 2.3 as a typical value. 
Then, one can plot the observed data in the left panels of Fig- 
ures [2] and [3] Placing the observed 2-10 keV photon index of 
3.06 ± 0.05 and the flux ratio R 2 ~ 6.8 ± 1.6 in the left panel 
of Fig. |2 we find that the data is marginally consistent with the- 
oretical prediction for the synchrotron radiation. However, the 
observed data point of 2-10 keV and 10-50 keV photon indices 
does not lie on any lines in Fig. [3] Furthermore, the observed 
data shows that for given value of 2-10 keV photon index, the 
10-50 keV photon index is smaller than expected via the syn- 
chrotron radiation. Hence one can claim the existence of an extra 
component. Note that, however, this argument is based on one- 
zone model. At present, it is uncertain whether 2-10 keV and 
10-80 keV X-rays come from the same emission region, which 
is left for future observations. 

One may insist on explaining the observed photon indices 
only by the synchrotron radiation. Indeed, it is marginally pos- 
sible that the observed 2-10 keV and 10-50 keV photon indices 
lie on the theoretical line of p > 3 and a ss 0.5. In this case, the 
acceleration is not limited by synchrotron cooling because gener- 
ally B > is expected (see Eq. [3J. This is apparently inconsistent 
with the observational implications that magnetic field strength 
is large enough (B > 0.5 mG) for the electron acceleration to be 
cooling limited dVink & Laming 200"l lBamba et alJl2005ah . 
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Appendix A: Characteristic energies 

In this section, we assume the Bohm diffusion, K(E) oc E, 
for sim plicity. More detailed analysis is found in Ohira et all 
d2012bl) . The extension to more general case {K(E) oc EP) is 
easy and omitted here. 
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Appendix A.1: Maximum energy of accelerated electrons 
and protons 

Let v s and be the shock velocity and the downstream mag- 
netic field, respectively. First, suppose that the maximum en- 
ergy of electrons is determined from the balance of the syn- 
chrotron loss and acceleration. Then, equating the acceleration 
time t- dcc (E) = IQ^cE /3ev 2 .Bd with the synchrotron cooling time, 
t syB (E) = 125 yr(£/10 TeV)- 1 (B d /100 ^G)" 2 , we obtain 



p(cool) 



24 



(to^)(i^g) " 2 TeV • (A - 1} 



where £ is a gyro factor. On the other hand, if the cooling is 
not significant, the maximum energy can be determined by the 
condition t icc (E) = f age , and we obtain 



7( a S e ) 



4.8 x 10 2 



( /' \ 2 (^-)(J^\ TeV . (A.2) 

£ UO^ms- 1 / yiOyuG/yiov/ 



For young SNRs, we typically expect the maximum energies of 
electrons and protons, £ malje and E maXtP are E^°°^ and E^lx, 



respectively. Then, using equation (IA. It and (IA.21 I. we eliminate 
v 2 /£ and obtain 



83 



10 TeV/ \l00fiG) \10 3 yr/ 



(A.3) 



Appendix A.2: The cooling break energy 



Let t be the characteristic time of SNR evolution, which may 
be the expansion time of f age itself. Then the cooling break, Eb, 
appears in the electron spectrum at the energy where f syn (£'b) = U 
that is, 



E b = 12.5 TeV 



t r 1 / B d - 



10 2 yr/ \100//G 



(A.4) 



The cooling break appears if Eb < E max<e = E^° l \ which can be 
rewritten as 

Bi>U9 ^^(w±^) (tov) (A - 5) 



Appendix B: Analytical derivation of asymptotic 
spectral shape 

Appendix B.1: The case of cooling-limited acceleration 

Zirakashv ili & Aharonianl (I2007I) obtained the asymptotic elec- 
tron spectrum near the maximum electron energy in the cooling- 
dominated phase, such as 



N(E)=A (E)exp[S (E)] 



(B.l) 



A (E) = £- 1/2 exp 



dE' 
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(B.2) 

Su(£) = -(£) 2 j- f (^*|^*) 2 ^, (B.3) 

where y s = 3r/(r - 1) (r is the shock compression ratio), and v s 
is the shock velocity. Functions K(E) and b(E) = dE/dt are the 
diffusion coefficient and energy loss rate of electrons, respec- 
tively. Subscripts 1 and 2 indicate upstream and downstream 
regions of the shock, respectively. In the case of synchrotron 
cooling, b(E) is proportional to E 2 . 

We assume that K(E) = KqEI 3 and b(E) = b E 2 where K 
and bo are constants, and that ratios b\/b2 and Kx/Kj are also 
constant. Then we obtain 

A Q (E)ccE l/2 . (B.4) 

However, for typical parameters of young SNRs, this positive 
slope can hardly be seen in the X-ray emission. Furthermore, ne- 
glecting non-dimensional terms on the order of unity, we derive 
Sq(E) sb -{Kobo/vl)E^ +l . By the way, acceleration time and 
synchrotron cooling time are calculated as t acc {E) « KoE^/v 2 
and f sy n(£) ~ Q)qE) , respectively, where we again neglect 
terms on the order of unity. Then equating them, we obtain 
^max 1 ' ~ (Vs/^o) 1/03+1) - Therefore, we finally derive 



So(E) 



y?+i 



p(cool) 
^max 



(B.5) 



Appendix B.2: The case of escape-limited acceleration 

We assume the test-particle regime and place a free escape 
boundary upstream at a distance of ( away from the shock front 
x = 0, that is, the particle distribution function is zero at x = 
—€ (< 0). Then stationary transport equation is s olved to find a 
particle spectrum around the shock front given by (Caprioli et al. 
l2009HReville et alJl2009h 

f 3r r E dlog£" 1 

TOcxexp , (B.6) 

{ r-l J 1 - exp [-v s {/K(E')] J 

where v s is the shock velocity. In the escape-limite d case, the 
maxim um energy is determined by the condition (Ohir aetaD 

l2010h 



K(E^) 



I . (B.7) 



Then, one can see K(E)/v s ( = K(E)/K(E ( ^) = (E/E^f. 
Changing variables into y - (E' /E^lt! we derive 

^{-^rm ■ 

where y(E) — {E/E i ^'f. In the case of E <k E ( ^ x \ the term 
e~ l ' y can be neglected, so that we obtain N(E) oc E- 3rl( - r ~ l) . On 
the other hand, if E » E^m, we approximate 1 - e~'/- v ~ l/y, 
resulting in N(E) oc exp[-(E / E^lx V 3 ]* where a term on the order 
of unity is again neglected. 



